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Model off-lattice sequences in two dimensions are constructed so that their native states are close to 
an on-lattice target. The Hamiltonian involves the Lennard- Jones and harmonic interactions. The 
native states of these sequences are determined with a high degree of certainty through Monte Carlo 
processes. The sequences are characterized thermodynamically and kinetically. It is shown that 
the rank-ordering-based scheme of the assignment of contact energies typically fails in off-lattice 
models even though it generates high stability of on-lattice sequences. Similar to the on-lattice 
case, Go-like modeling, in which the interaction potentials are restricted to the native contacts in 
a target shape, gives rise to good folding properties. Involving other contacts deteriorates these 
properties. 
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I. INTRODUCTION 

Understanding of the folding process of proteins 
is one of the current challenges in molecular bio- 
physics. Many insights into the nature of folding 
have been provided by studying lattice models in 
which a protein is represented by a chain of beads 
on a hupercubic lattice (see for instance ref. [Q). 
A more realistic modeling of proteins, however, re- 
quires considering off-lattice systems. Simple off- 
lattice systems have been discussed recently by Ir- 
back et al. and Klimov and Thirumalai The 
former authors have studied a model with two kinds 
of residues and they have found that very few such 
sequences give rise to a rapid folding. 

In this paper we focus on ways to design off-lattice 
sequences that form good folders. Specifically, we 
consider a two-dimensional (2D) target shape and 
assign two models of interaction energies between 
the beads. The first model is a generalization of the 
Go- like approach [|| to off-lattice situations. The 
second model, on the other hand, is a generaliza- 
tion of the rank-ordering-based assignment of the 
couplings 1^,1). When the beads are constrained to 
be located on sites of a lattice, both models pro- 
vide sequences which are good folders. Here, we 
demonstrate that this may not be true for off-lattice 
models. Namely, an extended character of the in- 
teractions, that is necessarily present in off-lattice 
Hamiltonians, gives rise to differing levels of frustra- 
tion in the two models and leaves only the Go-like 
sequences as good folders. This is compatible with 
the principle of minimum frustration proposed by 
Bryngelson and Wolynes 0|. 



We formulate our models in the context of the 
Hamiltonian used by lori, Marinari, and Parisi |^ 
(IMP) in which monomers interact via the Lennard- 
Jones potential. The amplitudes of the attrac- 
tive part between beads i and j of this poten- 
tial, Aij, representing a residue-dependent interac- 
tion, are quenched random variables. Additionally, 
the monomers are tethered sequentially along the 
chain by means of harmonic interactions. IMP and 
Struglia |^ have demonstrated the existence of a 
compact phase in a 3D version of the model. In these 
studies the dynamics have been defined in terms of 
a Monte Garlo process. The true ground state - the 
native state- however, typically is not known (as in 
IMP) or there is a substantial degree of uncertainty 
whether the state assumed to be native is indeed the 
ground state (as in the paper by Strugha ||]). 

Studies of the dynamics of folding require knowing 
the precise shape of the native state conformation 
without which the folding time cannot be defined. 
For small scale lattice models, the ground state may 
be obtained by an exact enumeration of conforma- 
tions. This method, however, does not apply to off- 
lattice models. Here we present 2D Lennard- Jones 
sequences of 16 monomers in which the ground state 
is known accurately and with a high degree of cer- 
tainty. We then provide some basic characteriza- 
tion of these sequences, obtained through the Monte 
Carlo procedure. 

The construction of the model sequences is pre- 
sented in Sec. II. We consider the lattice target 
shown in Fig. 1(a). In order to generate a Go- 
like sequence, G, we assign the Aij of the Lennard- 
Jones potential to be uniform in the native contacts 



and zero in the non-native contacts, i.e., all pairs 
of monomers that do not form contacts in the na- 
tive state interact only through a short-range repul- 
sion. The second model sequence, i?', is constructed 
by generating Gaussian Aij and adopting the rank- 
ordering scheme introduced recently ||^,^ in the con- 
text of lattice models. Sequence R' is an off-lattice 
analog of the sequence R discussed in Rcfs. [^,||. 
The principle here is to allocate the most strongly 
attractive Aij to those pairs of monomers which form 
contacts in the target compact lattice conformation. 
The true ground state of the sequence will necessar- 
ily be off-lattice but its shape will be close to the 
lattice target - it can be viewed as a somewhat dis- 
torted lattice target conformation. 

In Sec. Ill we describe the procedure adopted to 
determine the ground state. Basically, we form self- 
avoiding walks (SAW) in continuum space around all 
of the 69 maximally compact lattice conformations 
as the starting configuration of the sequence. We 
shall denote such starting SAW as CSAW to em- 
phasize closeness to compact structures. We then 
perform Monte Carlo quenches to reach local energy 
minima. The rank-ordered allocation of the cou- 
plings ensures that the lowest energy minima will 
be compact. For sequence R' we have shown that 
the lowest energy state is obtained from CSAW's 
which are near the target structure. We demon- 
strate that the use of arbitrary SAW's, i.e., which 
are not CSAW's, leads to local minima of higher en- 
ergies. For sequence G both the SAW's and CSAW's 
starting configurations easily lead to the native state 
which has the target compact shape. Our conclusion 
is that it is easier to find the native state for a good 
folder (G) compared to a bad folder (i?'). We have 
found also that only sufhciently large values of the 
spring constant k may guarantee the self-avoidance 
properties of the chain. 

We characterize geometries of conformations in 
terms of a certain distance away from a reference 
conformation. In Sec. IV we determine the proba- 
bility to stay in the basin and the folding tempera- 
ture. Both quantities are calculated by adopting a 
well defined characteristic basin size as obtained by 
the shape distortion approach In Sec. V we 

present results on the specific heat and structural 
susceptibility for R' and G and demonstrate that G 
is a good folder whereas R' is a bad folder. The 
bad folding properties of R' are in sharp contrast 
to what was found for the lattice sequence R 
and are due to the presence of many long-ranged 
couplings which impose conflicting constraints. Se- 
quence G is a good folder because there are many 



fewer constraints to satisfy. 

II. MODELS OF INTERACTIONS 

Following IMP we consider a self-interacting 
heteropolymer in 2D described by the Hamiltonian 
given by 

H = 5]{fc(d,., - do)%,j+i + 4[-^ ~ ^]} , 

(1) 

where i and j range from 1 to the number of beads, 
N, which in our model is equal to 16. The distance 
between the beads, dij is defined as \ri — rj\, where 
fi denotes the position of bead i. The harmonic 
term in the Hamiltonian with the spring constant 
fc, couples the adjacent beads along the chain. The 
remaining terms represent the Lennard- Jones poten- 
tial. In Isl Aij is chosen as Aij = Aq + y/erjij , where 
Aq is constant and rjij 's are Gaussian variables with 
zero mean and unit variance; e controls the strength 
of the quenched disorder. The case of rjij = and 
Aq — G would correspond to a homopolymer with 
the standard Lennard- Jones interaction used in sim- 
ulations of liquids. We adopt the units in which 
G=l and consider k to be either equal to 1 or to 25. 
We have found that the first choice, which has been 
used by IMP ||^, may lead to local energy minima 
in which the polymer is self-intersecting. 




(a) (b) 

FIG. 1. (a) Assigment of nine strongest couplings 
Aij to native contacts for sequence R' . The numbers 
indicate the relative strengths of the contacts, (b) The 
native conformation of sequence R' . 

For A^=16, there are 120 Lennard- Jones possible 
couplings between the monomers. The basic choice 
for the values of A^ that we shall use here is shown 



2 



in Tabic I. Here, all of the Aij's are positive, which 
corresponds to attraction. Table 1 also indicates a 
rank of a given attractive coupling if these are rank 
ordered from the strongest to the weakest attrac- 
tion. In nine of the couplings, the attraction is en- 
hanced by making the corresponding Aij bigger than 
one. These were chosen to coincide with the contacts 
present in the lattice target i?' shown in Fig. 1(a). 
The strongest attraction was assigned to be between 
beads 1 and 12 and the relative strengths of other 
native attractions are indicated in the figure. The 
remaining 111 couplings are assigned Aij with values 
which are smaller than 1. Overall the mean value Aq 
is equal to 0.784 and the dispersion to 0.205. 

The parameter do corresponding to the equilib- 
rium distance in the harmonic potential is chosen 
to be equal to 1.16, which is close to the equilib- 
rium position of the average Lennard-Jones poten- 
tial, {2C/Aoy^^. The target lattice shapes are built 
on a lattice with this lattice constant. For the Go- 
like sequence (G) we set Aij = 1 for native con- 
tact and for non-native ones. For this sequence 
we choose do = 2^1^ . The qualitative results do not 
depend, however, on the choice of do. 



III. LOCAL MINIMA AND NATIVE STATES 

In order to find spectra of the local energy minima 
we use the Monte Carlo procedure with local updat- 
ing moves. Monomers are moved randomly within 
a circle (the radius of circle varies from 0.0025(io to 
O.Oldo) away from their previous positions. Random 
quenching by starting from a SAW or a CSAW typ- 
ically yields a local energy minimum within of order 
2^* steps. In the minimum, the force acting on any 
of the monomers is at most of order 10~^ in units 
of the characteristic coupling (of order 1) divided by 
2C I Ai^)^!^ . Simulated annealing runs of comparable 
number of steps yielded similar results. The results 
reported on in this paper are based on the quenching 
procedure. 

The starting conformations were obtained by plac- 
ing monomers within a circle of radius 0.3c?q away 
from a SAW generated on the lattice or away from 
a maximally compact conformation on the lattice - 
CSAW. For each of the 69 maximally conformations, 
50 CSAW's were generated and the results were com- 
pared with those obtained based on 500 SAW's. 

We observe that there is a substantial gap, of or- 
der 4, between the minima obtained from CSAW's 
around the target and those obtained from all other 
starting configurations. 
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FIG. 2. Density of local minima for sequence B! 
and fc = 25. The left part corresponds to those ob- 
tained from deformed compact cells as starting confor- 
mations whereas the right part corresponds to the case 
when starting conformations are arbitrary self-avoiding 
chains. The energy gap is equal to 0.15. (For k = \ 
above E ~ —23.93 the distance between the beads 
may become bigger than the equilibrium distance of the 
Lennard-Jones potential by 10%. This happens in about 
1.4% of all conformations.) 

As to the choice of the elastic constant k we 
have found that fc=l is not very physical because 
the distances between consecutive monomers are not 
kept sufficiently rigid. Furthermore, we have ob- 
served that in states which are not local energy min- 
ima, the polymer conformations may become self- 
intersecting. This phenomenon was also observed in 
Ref. ||ll[ for shorter chains. Thus a stronger k is 
needed and we focus on fc=25. 

The energy histograms for the i?' and G sequences 
for k=25 are shown in Figs. 2 and 3, respectively. 
In the case of R' , the SAW configurations happen 
to yield results comparable to the CSAW configu- 
rations, but the statistical frequency of success in 
finding a low energy state favors the CSAW-based 
approach. The maximal distance between monomers 
does not exceed do by more than 4% in any local en- 
ergy minimum, which demonstrates very good self- 
avoiding properties. 

In the case of sequence G, about 25% of CSAW's 
trajectories and 12% of SAW's trajectories lead to 
the native conformation. Thus it is easier to find 
the native state for G than for R'. 
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The native state conformation of sequence R' is 
shown in Fig 1(b). The x ^ y coordinates of the 
beads of this conformation are hsted in Table 11. 
They are used in the studies of kinetics and thermo- 
dynamics properties. The native state of sequence G 
has exactly the lattice shape as shown in Fig. 1(a) 
because it is only the native contact energies here 
that one needs to minimize. 
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FIG. 3. The same as in Fig. 2 but for sequence G. 
The gap is about 2.9. 

As we shall show later, several tasks, like the de- 
termination of the folding temperature, are facili- 
tated by introducing a notion of a geometrical dis- 
tance between two conformations. A convenient def- 
inition of the distance 5ab between two conforma- 
tions a and 6 is provided by 



N 



(2) 



where f^' denotes the position of monomer i in con- 
formation a(h). The minimization is performed over 
translations, rotations and reflections. In practice, 
we put chain a over chain h by overlapping the two 
centers of mass, and then we find the optimal rota- 
tion of h which minimizes 5 ah- We pick the optimal 
angle from 1000 discrete values into which the 360° 
angle may be divided. 

The distances between the native conformations 
and their lattice counter parts, as in Fig. 1, are 
found to be equal to 0.296 and ^ for sequence R' 
and G respectively. 



IV. THE FOLDING TEMPERATURE 

We now proceed to the equilibrium characteriza- 
tion of the sequences introduced in this paper. The 
parameter that plays a primary role in determining 
the folding characteristics is the folding temperature 
Tf. In the case of lattice models, Tf may be defined 
m as the temperature at which the probability of 
occupying the native state becomes equal to \. For 
the off-lattice model, however, the native state has, 
strictly speaking, zero measure and one should deal 
with the probability of occupying the native valley. 
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FIG. 4. The dependence of Vo on T for sequence G 
and R' . The arrow indicates the position of the folding 
temperature which is equal to Tj — 0.23 for sequence G 
and Tf « 0.12 for R' . The results are averaged over 10 - 
40 Monte Carlo trajectories. 

One may define the folding temperature Tf 
through the following procedure. Suppose that S 
is the distance to the native state and 7-'((5) is the 
probability for a conformation to be in this distance 
away from the native state. Thus the probability 
to find the system in the immediate vicinity of the 
native state is given by 



Vo = 



V{5)d5 



(3) 
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where Sc is a cutoff distance. The folding tempera- 
ture Tf is then defined as the temperature at which 
^0 = 1/2. 

The size of the native basin Sc was estimated by 
the shape distortion approach [ pO| . In this approach 
one starts from the native state and performs ran- 
dom displacements of invidual beads in the chain, 
through a Monte Carlo routine. The distance to 
the native state is calculated by using Eq.(2) and 
the results are averaged over many Monte carlo tra- 
jectories. Below some critical temperature which 
may be interpreted as a folding temperature Tf, the 
distance to the native state gets saturated at suf- 
ficiently long time scales. The satuation value of 
this distance at the critical temperature can serve 
as the size of the native basin Sc- It is found that 
Sc « 0.2 and Sc ~ 0.09 for G and R', respectively 
pO| . The corresponding values of the folding tem- 
peratures obtained in this way are equal to Tf w 0.19 
and Tf « 0.09 for G and R' 0. 
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FIG. 5. The probability distribution function 7^(5) to 
find a state which is in distance S away from the native 
state for sequence R' . The values of T axe indicated next 
to the curves. 

Figure 4 shows the probability to be in the native 
valley versus temperature for G and R'. Vo has been 
determined by estimating the probability to stay in 
the native valley by starting in the native state and 
monitoring the system for 5x10® Monte Carlo steps 
(the distance updating is of order 0.01). Doubling 
the length of the run does not affect the results in 
any visible manner. For sequence G and R' , we have 



found Tf w 0.23 and Tf « 0.12. These values of 
Tf are close to those found by the shape distortion 
method @. 

Figure 5 shows the distribution 'P{S) at different 
temperatures for sequence R' . As the temperature 
increases the maximum becomes wider and moves 
toward larger distances to the native state S. 

We now turn to the temperature dependence of 
the folding time, tfoid- At high temperatures, reach- 
ing the native state takes long due to entropic effects. 
At low temperatures, on the other hand, glassy phe- 
nomena may set in and make the folding process 
extremely slow. Thus tfoid plotted against T typ- 
ically shows a minimum at a certain temperature 
Tmin- The idea of the existence of the glassy phase 
in proteins has originated in the Bryngelson and 
Wolynes studies of the random energy model (tHi^ 
and was subsequently tested in numerical simula- 
tions of lattice models by Socci and Onuchic [ p^ . 
Notice that Tf is a characteristic temperature that 
relates to equilibrium whereas Tmin is a character- 
istic temperature that relates to dynamical proper- 
ties. Even though at T = Tmin folding is the fastest 
this temperature also marks the onset of the glassy 
effects because they become stronger and stronger 
on departing from Tmin towards lower temperatures. 
Thus if Tf is significantly less than Tmin a sequences 
is bad folder. If Tf is comparable to Tmim or prefer- 
ably larger than Tmin, then the sequence is a good 
folder |1^. It should be noted that as a character- 
istic temperature that relates to dynamics one often 
uses the glass transition temperature Tg ||l|,^. Tg 
is, however, depends on a cutoff value of the time 
used in calculations. Our preference here is to use 
Tmim instead of Tg, not only because its definition is 
unique and independent of the value of the cutoff but 
also because the two quantities contain the same the 
physics: good folders are sequences in which glassy 
effect are not important around the folding temper- 
ature. 

In order to obtain tfoid we start from random 
configurations and evolve them through the Monte 
Carlo process until the native state is reached, tfoid 
is defined as the median number of MC steps af- 
ter which the system reaches the basin of the native 
state for the first time. The cutoff value of MC steps 
is taken to be 20 millions. The number of start- 
ing configurations varies from 20 to 40 depending 
on T. Our results are presented in Fig. 6 for se- 
quence R' and G. We have Tmin = 0.4 ±0.1 and 
Tmin = 0.15 ± 0.05 for R' and G respectively. Since 
Tf = 0.12 and Tf = 0.23 for these sequences one can 
see that R' is bad folder whereas G is a good one. 
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We emphasize that on-lattice Go and rank-ordering 
schemes lead to comparable ratios of Tf/Tmin- 
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FIG. 6. The temperature dependence of folding time 
for sequence G and R' . The results are obtained using 
20 - 50 starting configurations. 



V. THERMODYNAMIC PROPERTIES 

We now proceed to a further characterization of 
the sequences by considering the thermodynamic 
quantities such as the specific heat, C, and the struc- 
tural susceptibility, Xi tUl- It has been suggested 
^"^-|l6| that a small temperature difference between 
the maximum in C and the maximum in x when 
plotted against T indicates good folding properties. 
Thus studies of these two quantities may serve as a 
substitute for the information about the kinetics of 
folding. 

In the case of of an off-lattice model the departures 
of the sequence geometry from its native conforma- 
tion may be desribed through the structural overlap 
function [S| as: 



X. = 1- 



1 



(4) 



where rij is the distance between the beads i and 
j for a given conformation, is the corresponding 
distance in the native conformation, and Q{x) is the 
Heavyside function. Here 6c denotes the size of basin 
as defined above in the previous Section. If |r.y — 
^ <5c the beads i and j are assumed to form a 
contact within the native valley. The susceptibility- 
like parameter x is defined as the thermal fluctuation 
of Xs- 



X{T) - < xl{T) > - < Xs(T) >2 



(5) 



where the angular brackets indicate a thermody- 
namic average. 
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FIG. 7. The temperature dependence of x S'ld C for 
G and R' . The starting configurations are the native 
state. The results are averaged over 10 MC trajectories. 
The arrows indicate the position of T/ found from the 
condition that the probability to stay in the native state 
is equal to 1/2. 

The specific heat is defined in the usual way, i.e. 
by the energy fluctuations: 
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A peak in C may be interpreted as corresponding 
to the onset of slow kinetics and the peak in x 
corresponding to the folding temperature. 
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We calculate the thermodynamic quantities using 
the native state as the starting configurations. The 
results are averaged over many MC trajectories. The 
equilibration is checked by monitoring the stability 
of the data against at least three-times longer runs. 
We have used 5x 10^ MC steps as in the studies of Vo 
(the first 2.5x10^ steps are not taken into account 
when averaging). This number of MCS appears to 
be enough to equilibrate the system in the temper- 
ature interval we use. We have also found that it 
is harder to reach the equilibrium using SAW's as 
starting configurations. 

The temperature dependence of x s^nd C for se- 
quence G and R' is shown in Fig. 7. The peak in 
X almost coincides with Tf. For G maxima of x 
and C are located at the same position suggesting 
that sequence G is a good folder Q. This results 
agrees with that obtained in the previous section by 
studying kinetics of folding. 

For sequence R' , the peak in C is found to be 
broad and it is located at a substantially higher T 
than the peak in x confirming bad folding properties 
of R' . So studies of thermodynamics properties and 
folding kinetics lead to the same conclusion about 
folding in our model. 



VI. CONCLUSIONS 

In this work, we have constructed and studied 
two off-lattice Lennard- Jones models in two dimen- 
sions. The long-range nature of interactions between 
amino-acids in off-lattice models may lead to too 
much frustration, as in the case of sequence R' , and 
thus to bad foldicities. Restricting the number of in- 
teractions only to the native contacts, as in sequence 
G, may reduce the frustration and bring about good 
folding properties. We have demonstrated this by 
studying both thermodynamic and dynamical prop- 
erties of the two sequences. The values of the folding 
temperature found from the probability to get out 
from the native state roughly agree with those ob- 
tained by the shape distortion approach. We have 
studied the kinetic of folding with a simplified MC 
dynamics whereas previous works have been focused 
on equilibrium aspects. Similar to the on-lattice 
models, both approaches give the same information 
about folding properties. 

We thank J. R. Banavar for discussions. This 
work was supported by Komitet Badan Naukowych 
(Poland; Grant number 2P03B-2513). 
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10 


0.802427 


48 


3 


13 


0.518711 


116 


6 


10 


0.676479 


76 


9 


16 


0.850939 


39 


1 


11 


0.763514 


55 


3 


14 


0.660732 


85 


6 


11 


0.567934 


109 


10 


11 


0.935812 


22 


1 


12 


1.498411 


1 


3 


15 


0.852925 


38 


6 


12 


0.571524 


108 


10 


12 


0.951547 


16 


1 


13 


0.899355 


30 


3 


16 


1.050544 


9 


6 


13 


0.658820 


86 


10 


13 


0.521468 


114 


1 


14 


1.396376 


3 


4 


5 


0.878341 


33 


6 


14 


0.512926 


119 


10 


14 


0.930277 


25 


1 


15 


0.723510 


66 


4 


6 


0.763190 


56 


6 


15 


0.612393 


99 


10 


15 


0.889378 


31 


1 


16 


0.655780 


87 


4 


7 


1.269149 


8 


6 


16 


0.628334 


95 


10 


16 


0.590144 


104 


2 


3 


0.720592 


68 


4 


8 


0.577229 


106 


7 


8 


0.670401 


78 


11 


12 


0.932459 


23 


2 


4 


0.541154 


112 


4 


9 


0.614099 


96 


7 


9 


0.647826 


89 


11 


13 


0.528390 


113 


2 


5 


0.770339 


53 


4 


10 


0.918839 


27 


7 


10 


0.642560 


92 


11 


14 


0.921718 


26 


2 


6 


0.938498 


21 


4 


11 


0.597736 


102 


7 


11 


0.965115 


12 


11 


15 


0.770826 


52 


2 


7 


1.291634 


7 


4 


12 


0.571611 


107 


7 


12 


0.834125 


42 


11 


16 


0.955652 


14 


2 


8 


0.696944 


72 


4 


13 


0.561490 


110 


7 


13 


0.753538 


58 


12 


13 


0.629158 


94 


2 


9 


0.603317 


100 


4 


14 


0.725428 


64 


7 


14 


0.561397 


111 


12 


14 


0.884585 


32 


2 


10 


0.832387 


43 


4 


15 


0.954825 


15 


7 


15 


0.712665 


70 


12 


15 


0.949139 


18 


2 


11 


0.672241 


77 


4 


16 


0.511493 


120 


7 


16 


0.864981 


35 


12 


16 


0.932023 


24 


2 


12 


0.821738 


46 


5 


6 


0.687186 


74 


8 


9 


0.943903 


20 


13 


14 


0.848886 


40 


2 


13 


0.644946 


90 


5 


7 


0.665046 


81 


8 


10 


0.660797 


84 


13 


15 


0.613872 


97 


2 


14 


0.689698 


73 


5 


8 


0.686304 


75 


8 


11 


1.301709 


6 


13 


16 


0.734204 


63 


2 


15 


1.409250 


2 


5 


9 


0.771102 


51 


8 


12 


0.903380 


29 


14 


15 


0.519439 


115 


2 


16 


0.655129 


88 


5 


10 


0.595747 


103 


8 


13 


0.612414 


98 


14 


16 


0.963250 


13 


3 


4 


0.781090 


49 


5 


11 


0.707574 


71 


8 


14 


0.825123 


44 


15 


16 


0.866528 


34 



TABLE I. Values of the Aij for sequence R' . The rank of a coupUng is indicated to the right of its value. 
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monomer 


Xn 


VN 


monomer 


Xn 


VN 


1 








9 


-0.037572 


2.225111 


2 


1.005509 


0.573347 


10 


-1.073118 


2.718455 


3 


2.091357 


0.126744 


11 


-0.992672 


1.573907 


4 


2.009697 


1.279681 


12 


-0.956188 


0.436143 


5 


1.946778 


2.424499 


13 


-0.972265 


-0.723659 


6 


0.845883 


2.806329 


14 


0.116648 


-1.062284 


7 


0.997773 


1.647486 


15 


1.123476 


-0.470181 


8 


-0.035144 


1.079819 


16 


2.154153 


-0.984088 



TABLE II. The x — y coordinates of the native conformation for R' and k = 25. 
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